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Abstract 


A  computer  system  serving  both  batch  and  interactive  jobs  is  modeled 
as  a  single  server  queue,  with  an  infinite  source  of  ordinary  customers  and 
a  finite  source  of  priority  customers. 

The  stability  condition  and  the  stationary  probability  distribution 
are  determined.  For  a  stable  system,  it  is  shown  how  the  distribution  and 
moments,  of  the  waiting  time,  the  sojourn  time,  the  completion  time,  and 
the  busy  period,  may  be  efficiently  computed. 


Keywords 

Queueing  system,  finite  priority  source,  batch  and  interactive 
computer  model,  computational  probability. 
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Introduction 

The  purpose  of  this  paper  is  to  present  the  algorithmic  analysis  of 
a  single  model  for  a  computer  system,  serving  batch  and  interactive  jobs. 

The  methodology  used  here  has  been  introduced  in  recent  years  by  Neuts  [8-12]. 

The  computer  system  under  study  is  represented  by  a  single  server 
queue  with  infinite  waiting  room  (see  Figure  1).  Batch  jobs  arrive  according 
to  a  Poisson  process  with  parameter  x^.  The  time  needed  to  process  a  batch 
job  is  exponential  with  parameter  Once  a  batch  job  is  terminated,  it 
leaves  the  system.  Interactive  jobs  are  submitted  by  a  finite  number,  N, 
of  interactive  terminals.  The  time  needed  by  a  terminal  to  submit  a  job 
is  exponential  with  parameter  X2_  Once  the  job  is  submitted,  the  terminal 
must  wait  until  the  job  is  processed,  at  which  time  the  terminal  initiates 
the  submission  of  a  new  job.  The  time  needed  by  the  server  to  process  an 
interactive  job  is  exponential  with  parameter  The  interactive  jobs 
have  preemptive  priority  over  the  batch  jobs.  It  will  readily  be  observed 
that  such  a  system  is  a  special  case  of  the  priority  queueing  model,  with 
infinite  ordinary  source  and  finite  priority  source,  which  has  already  been 
studied  by  Avi-Itzhak  and  Naor  [1],  Jaiswal  and  Thiruvengadam  [3], 
Thiruvengadam  [13,14],  and  most  comprehensively  by  Jaiswal  [4].  The 
results  in  [4]  are  obtained  via  the  analysis  of  the  renewal  process  formed 
by  successive  busy  periods  of  the  queueing  system.  The  results  mainly 
consist  of  Laplace- Stieltjes  transforms  of  distributions  or  of  generating 
functions.  These  do  not  simplify  very  much  in  the  single  case  considered 
here,  where  the  service  time  distributions  are  exponential. 

In  contrast,  the  method  used  here  will  yield  computationally  efficient 
expressions.  We  shall  determine  the  stability  condition  for  the  system  and 
prove  that  the  stationary  probability  vector  has  a  matrix-geometric  form. 
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Appropriately  partitioning  that  vector  x  as  (Xq»x.|,...),  it  is  shown  that 

4,  -  l(I-R)R\  for  i*0, 

where  the  positive  vector  £,  and  the  positive  matrix  R  are  explicitly 
determined.  We  show  in  Sections  2  to  4  how  the  distribution  and  moments 
of  the  waiting  time,  the  sojourn  time,  the  completion  time,  and  the  busy 
period  are  determined.  A  comprehensive  discussion  of  stochastic  models 
with  embedded  Markov  chains  having  a  matrix-geometric  invariant  probability 
vector  may  be  found  in  Neuts  [12]. 

To  conclude  this  introduction,  we  remark  that  the  simple  network 
described  in  Figure  1  is  not  amenable  to  the  approach  of  Baskett  et.al.  [2], 
since  the  two  priority  classes  of  jobs  have  different  service  rates. 

1.  The  Stationary  Distribution 

Under  the  assumptions  of  Poisson  arrivals  and  exponential  services, 
the  model  may  be  described  as  a  continuous-parameter  Markov  chain  on  the 
state  space  {(i,j),  i^O,  Osj <N > ,  where  i  and  j  respectively  denote  the 
number  of  batch  and  interactive  jobs  in  the  queueing  system. 

The  infinitesimal  generator  Q  of  the  Markov  chain  is  a  block-tridiagonal 
matrix  of  the  form 


where  AQ,  A^  and  Ag  are  square  matrices  of  order  N+l ,  defined  as  follows. 
The  matrix  AQ  corresponds  to  transitions  from  states  (i,j)  to  (1+l,j')» 
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for  i^O,  and  is  equal  to  x^l.  The  matrix  A2  corresponds  to  transitions 
from  states  (i,j)  to  for  1^1 ;  (^)qq  is  equal  to  ur  all  the 

other  elements  of  A2  are  equal  to  zero.  The  matrix  Aj  corresponds  to 
transitions  from  states  (i,j)  to  (i,j'),  for  1^1 ,  and  Is  given  by 

-x1-NX2-m1  NX 2  0 

w2  -Xr(N-l)X2-u2  (N-l  )X2 

A  0  u2  -X^tN-ZjXg-Ug 


Let  denote  the  vector  of  stationary  probabilities.  It  is  the  unique 
solution  to  the  system  xQ=M),  >^e=l,  where  e  represents  a  vector  with  every 
component  equal  to  1.  Let  us  partition  x  as  (xg.x^ ,. . . ) ,  where  )(..  is  an 
N+l-vector,  and  corresponds  to  the  states  {(i,j),  0<j <N } .  Furthermore,  we 
define  the  matrix  A  by  A=Aq+A^+A2>  It  is  the  infinitesimal  generator  of  a 
finite,  irreducible,  continuous-parameter  Markov  chain.  By  *  we  denote 
its  vector  of  stationary  probabilities,  i.e.  *Aa0,  ^e=l. 

Lemma  1 

The  vector  is  given  by 

ni  i  [  N  mi 

"j =  mrr  °2y0  twit  p2j  ’  for  0sJiN’ 

where  p2=x2/u. 

Proof.  The  proof  is  elementary. 

Theorem  1 

The  system  is  stable  if  and  only  if 
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P  =  Xl  Ck^o  TN^T!"  p2^  /Wl  <  (2) 

If  the  system  is  stable,  the  stationary  probability  vector  x  is  given  by 
x=(xq,*j,...)»  where 

x.  =  *(I-R)r\  for  i^O.  (3) 

The  matrix  R  is  the  unique  nonnegative  solution  with  maximal  eigenvalue 
strictly  less  than  one  of  the  matrix  quadratic  equation 

R2A2+RA1+Aq  =0.  (4) 

The  matrix  R  is  equal  to  lim  Rn,  where  the  matrices  Rp  are  defined  as 

n->« 

follows. 


Rn+1  =  ■A0Al’1*RnA2Al’1’  for  n>0.  (5) 

Moreover,  the  matrix  R  is  strictly  positive,  and  its  first  column  is  equal 
to  p-je,  where  p-j=x^/y ^ . 

Proof.  This  theorem  is  proved  by  repeating  arguments  given  in  [8-11].  We 
refer  the  reader  to  those  papers  for  details  and  only  indicate  the  main 
steps  of  the  proof  here. 

Firstly,  it  results  from  [11],  Theorem  1,  that  if  the  system  is  stable, 
then  x^XqR1*  for  i*0,  where  the  matrix  R  satisfies  the  equation  (4),  the 
maximal  eigenvalue  of  R  is  strictly  less  than  one,  the  matrix  R  is  the 
minimal  nonnegative  solution  of  (4). 

Secondly,  we  show  that  the  sequence  of  matrices  Rn  defined  in  (5), 
converges  monotonically  to  the  minimal  nonnegative  solution  of  (4).  Also, 
we  show  that  Rn  is  strictly  positive,  for  nsl. 


Thirdly,  we  prove  that  the  matrix  quadratic  equation  (4)  has  a 
unique  solution  with  maximal  eigenvalue  strictly  less  than  one  if  and  only 


if  JL^e^Age.  It  results  from  Lemma  1  that  the  stability  condition 
n/^e^Age  is  equivalent  to  the  inequality  (2).  It  is  easy  to  verify  that 
if  the  system  is  stable,  then  Xg=]^(I-R),  which  completes  the  proof  of  (3). 
Finally,  we  show  that  RA^Age,  which  proves  that  the  first  column  of  R  is 
equal  to  p-je. 

In  general,  for  systems  which  may  be  analyzed  by  the  approach  followed 
here,  the  matrix  R  must  be  numerically  computed  by  using  the  recurrence 
relation  (5).  In  our  case  however,  it  is  possible  to  go  further  in  the 
determination  of  the  matrix  R,  as  we  show  in  the  next  theorem. 

Theorem  2 

The  matrix  R  is  explicitly  given  by 

x2 

R  ■ -Vi'1  ♦  v'"-  (6) 

where  m  denotes  the  first  row  of  the  matrix  A,"\  and  the  matrix  M  is  equal 
to  e*m,  i.e.  each  row  of  the  matrix  M  is  equal  to  m. 

Proof.  The  matrix  R  satisfies  the  equation 

R  =  -aoai"1"R2a2ai"1.  (?) 

From  the  structure  of  Ag  and  it  results  that  (7)  may  be  written  as 

R  =  -x1A1“1-ii1R2H*,  (8) 

where  the  first  row  of  M*  is  equal  to  m,  the  remaining  rows  are  equal  to  0. 
Clearly,  (8)  is  equivalent  to 
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2 

where  r*  denotes  the  first  column  of  R  ,  the  first  column  of  R*  is  equal 
to  r*,  the  remaining  columns  are  equal  to  0.  Since  the  first  column  of  R 
is  equal  to  (by  Theorem  1),  then  r*=p^Re,  and  (9)  becomes 

R  =  -x1A1"1-X1Re*m.  (10) 

Postmulti plying  both  sides  of  (10)  by  e,  and  replacing  the  obtained  value 
for  Re  in  (10),  completes  the  proof  of  the  theorem. 

It  is  now  a  simple  matter  to  determine  special  probabilities  of 
interest.  Let  £  and  £  respectively  represent  the  marginal  distribution  of 
the  number  of  interactive  and  batch  jobs.  Clearly,  one  has  that 

co 

£  =  z  x,  =  1-  The  vector  w  is  explicitly  given  in  Lemma  1,  and  there  is 
i=0 

no  difficulty  in  computing  the  moments  for  the  marginal  distribution  of 
the  interactive  jobs.  To  determine  the  distribution  and  moments  for  the 
batch  jobs  is  only  slightly  more  difficult.  We  shall  only  state  the 
following  result.  The  proof  is  identical  to  that  of  [5],  Lemma  1,  and  is 
not  presented  here. 


Lemma  2 

The  marginal  distribution  of  the  number  of  batch  jobs  is  given  by 


p.j  =  ^(I-RjR’e,  for  i^O. 

The  moment  3V,  about  the  origin,  is  equal  to  B  =itX  e ,  for  v>l ,  where 
the  matrices  Xv  of  order  N+l  are  recursively  defined  as  follows. 


-1 


X0  =  I.  Xn+1  =  R(I-R)*  Z 


j=0 


for  nsO. 


In  particular,  one  has  that 
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@1  =  wR(I-R)_1e, 

62  =  iR(I+R)(I-R)-2e, 
and  e3  =  *R(I+/,R+R2)(I-R)‘3e. 

Finally,  we  determine  easily  the  following  probabilities  for  the  state 
of  the  server.  The  probability  that  the  server  is  idle  is  equal  to  Xg  o=no"pl* 
The  probability  that  the  server  is  busy  and  processing  a  batch  job  is  equal 

co 

to  z  x-  n=p,.  The  probability  that  the  server  is  busy  and  processing  an 
i=l  1,u  1 

interactive  job  is  equal  to  l-Pg=l-iTg. 

2.  The  Waiting  Time  and  the  Sojourn  Time  for  a  Batch  Job 

It  is  obvious  that  the  interactive  jobs  actually  form  an  M/M/1  queue 
with  finite  source.  This  type  of  queue  has  been  widely  studied  and  therefore, 
we  will  only  examine  the  batch  jobs.  This  can  be  done  by  any  of  two 
methods.  Firstly,  by  using  the  concept  of  phase  type  distributions  -  Neuts 
[6,7]  -  we  obtain  an  explicit  form  for  the  probability  distributions. 

The  second  approach  provides  expressions  which  are  better  suited  to 
numerical  computations. 

We  consider  the  queues  in  steady- state  and  we  denote  by  w  the  time 
spent  by  a  batch  job  in  the  incoming  queue,  by  Wq(x,i),  the  conditional 
distribution  of  w,  given  that  upon  arrival  of  the  job,  the  system  contained 
i  batch  jobs  already  and  by  W  (x)  the  distribution  of  w. 

M 

We  will  need  the  following  matrices  and  vectors:  Sg  is  the  square 
matrix  of  order  N  obtained  by  deleting  the  first  row  and  the  first  column 
of  A-j+Ag,  S.j  (i>0)  is  the  square  matrix  of  order  i(N+l)  given  by: 
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where  C-j*A^+Aq  and  Cq=A2;  is  the  N-vector  obtained  by  deleting  the 
first  component  of  x^,  and  a..  (i>0)  is  the  i  (N+l  )-vector  (0 ,0, . . .  ,0,)^ ) . 

Theorem  3 

The  distribution  of  w  is  given  by: 

oo 

W„(x)  =  1  -  z  a.  exp(S..x)e,  for  x>0. 

q  i=0  1  1 

The  moments  are  given  by: 

oo 

E[wV]  =  ( - 1 ) V •  v  1  z  a.S/Ve,  for  v>0. 
i=0  1  1 

Proof.  Given  that  the  system  contains  i  batch  jobs  already,  w  is  equal  to 
the  time  until  absorption  in  a  continuous  parameter  Markov  chain,  with 
transitions  among  transient  states  governed  by  ,  and  initial  probability 
vector  .  Therefore 

Wq(x,i)  =  l-tx^e)"1^  exp^xje,  for  x>0, 

and  E[wv|i]  =  (-l)v  vKx.^'VSi'V 

The  proof  is  now  immediate. 

Remark.  Although  this  theorem  provides  us  with  an  explicit  expression  for 

W  (x),  it  is  difficult  to  use  in  order  to  get  numerical  results,  because  it 

q  /  cn„n 

00  /  ®  S.x 

contains  a  double  series  z  «,(  z  ■■■  ■ 

i=0  1 \n=0  n- 

These  series  are  convergent,  but  involve  large  matrices.  Nevertheless, 
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the  distribution  of  w  may  numerically  be  obtained  by  another  method. 

Clearly,  w  is  the  time  until  absorption  in  the  Markov  chain  with 

infinitesimal  generator  Q  and  initial  probability  vector 

w 

xf(xo*-l*“**^k**“)*  where 


and  C  is  the  matrix  with  a  first  row  identically  zero. 

Let  y.  .(x)  be  the  probability  that  this  Markov  chain  is  in  the  state 

*  J 

(i,j)  at  time  x.  We  have  that  Wq(x)=yQQ(x).  The  Kolmogorov  equations  for 
this  chain  are 

^(x)  =  yQ(x)C+x1  (x)C0,  ^ 

Zj(x)  =  Zi(x)C1+yi+1(x)C0,  for  1*1, 

with  the  initial  conditions  y. (0)*x.j • 

The  Markov  chain  Qw  can  only  move  towards  lower  states.  It  is  there¬ 
fore  obvious  how  to  truncate  the  infinite  system  of  differential  equations 
(11).  In  order  to  lose  a  probability  mass  of  at  most  e  in  the  tail  of  the 
distribution  function  W^(x),  one  truncates  at  the  index  K  such  that 

CO 

I  ^(I-R)Rve  =  jtR^*  e  <  e 
v=K+l 

This  approach  does  not  lead  to  an  explicit  form  for  W^(x),  but  is 
easily  implemented. 

We  define  the  sojourn  time  of  a  batch  job  in  the  system  as  the  time 
between  the  moment  the  job  enters  the  system  and  the  moment  it  leaves  the 
system.  The  sojourn  time  can  be  analyzed  exactly  like  the  waiting  time. 
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3.  The  Completion  Time  of  a  Batch  Job 

Let  w  (x)  denote  the  distribution  function  of  the  completion  time  w_ 
of  a  batch  job,  i.e.  the  time  between  the  moment  when  the  processing  of 
the  job  begins,  and  the  moment  when  the  job  leaves  the  system.  wc  is  also 
the  time  until  absorption  into  the  absorbing  state  of  the  Markov  chain 
with  (N+2)  states  and  infinitesimal  generator  Qc  given  by 


0  10  ...  01 
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and  with  the  initial  probability  vector  v=(0,en)=(0,l ,0,. .. ,0).  Clearly 
now, 

wc(x)  =  1-e0  exp(C1x)e, 
and  E[wp  =  (-l)v  v!e0C1~ve, 

Remarks . 

a.  E[w  ]  is  equal  to  f-  =  . 

Indeed,  E[wc]  is  equal  to  minus  the  sum  of  the  elements  of  the  first 
row  of  C,”^;  these  elements  are  easily  computed,  since  they  are  the  solution 
of  a  system  that  differs  only  from  the  system  for  w  in  the  first  equation. 

It  is  now  easy  to  prove  (13). 

b.  Equation  (12)  may  be  written  as 

W  (x)  =  -  z  c.xk,  (14) 

c  k=l  K 


for  x>0, 
for  v>l . 


(12) 


(13) 
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1  k 

where  is  the  sum  of  the  elements  of  the  first  row  of  C-j.  It  is  not 

necessary  however  to  compute  all  the  components  of  C-j  in  order  to  obtain  c^. 
/  k)  Ik 

If  we  denote  by  c/  1  the  first  row  of  Cj1,  one  has  that 


c(k)  =  1  c(k-l)r 
-  k-  V 


and 


r  _  c(k)  .  (k-1 ) 

ck  '  -  -  "  '  k  ^  V 


c.  The  series  (14)  is  unfortunately  difficult  to  use  for  numerical 

purposes  because  its  terms  alternate  in  sign  and  have  large  absolute 

30 

values  (greater  than  10  in  some  cases).  Since  this  series  has  a  sum  in 
the  interval  (0,1),  it  is  clear  that  it  is  computationally  unstable. 
Nevertheless,  one  easily  obtains  the  distribution  of  wc  as  the  solution 
of  a  finite  system  of  differential  equations  by  using  the  same  approach  as 
for  the  waiting  time. 


4.  The  Busy  Period  of  the  System 

The  busy  period  is  defined  as  the  interval  of  time  between  the  moment 
when  the  CPU  becomes  active  (by  the  arrival  of  either  an  interactive  or  a 
batch  job)  and  the  first  moment  when  the  CPU  becomes  again  inactive.  We 
will  have  to  distinguish  two  different  types  of  busy  periods,  according  as 
it  begins  with  the  arrival  of  an  interactive  (type  1)  or  a  batch  job  (type  0). 

We  denote  by  hgfk.fc.x),  ksO,  £>0,  the  probability  that  starting  in  the 
state  (i+1,0),  i^O,  at  time  0,  in  the  Markov  chain  with  infinitesimal 
generator  Q,  the  first  visit  to  the  state  (i,0)  occurs  no  later  than  time  x, 
and  that  exactly  k  batch  and  i  interactive  jobs  are  processed  during  that 
first  passage  time.  Similarly,  we  define  h.(k,t,x),  j=l,2,...,N,  for  the 

J 

first  passage  time  from  the  state  (i,j),  i >0 ,  to  the  state  ( i ,0) .  Further¬ 
more,  we  introduce  the  transforms 
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h*(z,y,s)  =  z  E  zky*  /"  e'sxdh.(k,i,x). 
J  k=0  *=0  0  J 


H*(z,y,s)  is  a  square  matrix  of  order  (N+l)  such  that  HJ0(z,y,s)=hJ(z,y,s) 
for  j=0,l,...,N,  the  remaining  elements  are  all  zero  and  fi=H*(l ,1 ,0). 
Finally,  A£  is  a  square  matrix  of  order  (N+l)  with  only  one  non-zero  element 
(Apio  =  n2- 

By  following  the  argument  in  [9],  we  prove  that  the  matrix  H  is 
stochastic  if  the  system  is  stable.  Then  each  row  of  H  is  equal  to  e^. 

Let  us  denote  by  m,  ng,  the  following  vectors 


m  = 


'  Is  H*(z,y,s)ej 
L3s  J(z=l  ,y=l  ,s=0) 

:  f|r  H*(z,y,s)e1 
L3Z  J(z=l,y=l,s=0) 

:  f-|vH*(z,y,s)e| 

L9y  J(z=l ,y=l ,s=0) 


Their  components  mp,  nBp,  nIp,  for  p=0,l,  are  respectively  equal  to  the 

mean  busy  period  and  the  mean  numbers  of  batch  and  interactive  jobs 
processed  during  a  busy  period  of  type  p. 


Theorem  4 

If  the  queueing  system  is  stable  (p<1),  the  vectors  m,  rig,  and  jij 
are  the  solutions  of 


Mm  =  -e 

(15) 

i  sir 

i 

ii 

(16) 

=£1* 

i 

H 

# 

(17) 

where  M  is  a  non- singular  matrix  given  by 


M  -  A1+A0*A0H-A'  +  A2 
and  Ei  =  t(p-)  ,0,0,. . .  ,0) 


H-2  “  (0»vJ2*0».  • ,  *0) 


One  obtains  the  following  explicit  forms: 


m0  =  Pf0(l-P)  *  ml  =  "0NA2(1-p)  * 


'bo  =  TT^py  ’ 


.  M1^ 


B1 


nI0  "  w^l-P)  ’ 


Proof.  The  proof  of  this  theorem  is  purely  technical  and  may  be  found  in 


the  appendix. 


5.  Other  Models 


In  this  section  we  shall  briefly  examine  three  queueing  systems 
closely  related  to  the  system  defined  in  the  introduction. 

Firstly,  we  shall  consider  the  case  where  interactive  jobs  have 
non-preempti ve  priority  over  the  batch  jobs.  If  an  interactive  job  is 
submitted  while  a  batch  job  is  being  processed,  the  interactive  jobs  in  the 
queue  must  wait  until  that  job  leaves  the  system.  In  this  case,  the  state 
space  may  be  represented  by  ((i,j),  i=0,  Osj<N;  i>0,  -N^j<N},  where  j>0 
indicates  that  a  batch  job  is  being  processed,  while  j<0  indicates  that  a 
batch  job  is  being  processed,  and  |j|  interactive  jobs  are  waiting.  The 
infinitesimal  generator  Q  is  now  given  by 
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Q 


B0  0 

A1  A0 
A2  A1 


(21) 


where  the  square  matrices  Aq,  A-j  ,  and  ^  are  order  2N+1,  the  square 
matrix  is  of  order  N+l,  the  matrices  Bq  and  ^  are  rectangular,  with 
the  appropriate  dimensions.  Each  of  these  matrices  is  very  sparse,  we 
do  not  indicate  their  structure,  since  it  is  cumbersome  and  would  serve  no 
purpose  here.  It  is  clear  that  this  model  may  be  solved  by  using  exactly 
the  same  approach  as  in  Sections  1  through  4. 

Secondly,  we  may  consider  the  preemptive-del ayed  priority  discipline. 

If  an  interactive  job  is  submitted  while  a  batch  job  is  being  processed,  the 
latter  is  allocated  a  quantum  of  time,  say  D,  during  which  it  may  terminate 
its  processing.  If  it  does  not  leave  the  system  before  D,  it  is  sent  back 
in  the  waiting  queue  and  the  interactive  jobs  are  processed.  If  we  assume 
that  D  has  an  exponential  distribution,  then  the  system  is  very  similar  to 
the  non-preemptive  priority  models.  Only  the  matrix  A-j  in  (21)  is  modified. 

Thirdly,  let  us  assume  that  the  system  contains  c>l  servers.  Then 
some  modifications  occur,  and  we  briefly  present  the  most  important  ones. 

For  each  of  the  three  queueing  disciplines,  the  structure  of  the  matrix  Q 
is  slightly  modified.  In  addition,  for  the  preemptive-del ayed  and  the 
non-preemptive  disciplines,  a  new  variable  must  be  introduced  to  describe 
the  mixture  of  batch  and  interactive  jobs  being  processed.  This  adds  to 
the  dimensionality  of  the  problem  since  the  elements  themselves  of  the 
matrices  A  and  B  in  (21)  become  matrices.  For  the  preemptive  discipline, 
no  new  variable  is  necessary  and  the  stationary  distribution,  the  waiting 
time,  and  the  busy  period  may  be  analyzed  in  much  the  same  way  as  in 
Sections  1,  2,  and  4.  However,  the  analysis  of  the  completion  time  of  a 
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batch  job  will  be  different.  It  will  be  necessary  to  introduce  a  rule 
for  selecting  which  batch  job  is  to  be  interrupted.  Different  rules  will 
yield  different  distributions  for  the  completion  time.  It  is  our  intention 
to  discuss  this  point  elsewhere. 

6.  Numerical  Results 

The  parameters  used  in  this  section  have  been  obtained  after  the 
analysis  of  some  statistics  from  the  Computer  Center  at  our  University; 
we  obtained  the  following  values: 

A i  =  0.030  sec'1 
=  0.056  sec  1 
u-j  =  0.5  sec'1 

^2  =  2.4  sec-1 

N  =  43  (mean  number  of  active  terminals). 

We  used  four  of  these  values  and  let  the  fifth  one  vary.  We  will 
show  here  some  results  obtained  by  varying  the  number  of  terminals  or  the 
interactive  service  rate. 

We  show  in  Figure  2  how  the  traffic  coefficient  p (n )  depends  on  the 
number  of  terminals.  It  may  be  observed  that  43  active  terminals  is  a 
number  near  the  critical  region;  a  few  supplementary  terminals  will  bring 
the  system  to  saturation. 

Table  I  gives,  for  different  values  of  N,  the  following  information: 

-  the  traffic  coefficient  p  defined  in  Theorem  1; 

-  the  mean  number  E[j]  of  interactive  jobs  and  the  percentiles 
.5,  .9,  .99,  and  .999,  of  the  distribution  £; 

-  the  mean  number  E[i]_of  batch  jobs  and  the  same  percentiles 
of  the  distribution  j>; 


ar-f 


-  the  mean  value  E[Wj]  and  the  percentiles  .5,  .9,  .99  and  .999, 
for  the  waiting  time  distribution  of  interactive  jobs; 

-  the  expected  waiting  time  E[w]  and  completion  time  E[wc]  for 
batch  jobs,  the  mean  number  mjc  of  interactive  jobs  processed 
during  the  completion  time  of  a  batch  job; 

-  the  expected  length  mn  and  m,  of  a  busy  period  of  type  0  and 

type  1.  u  i 

On  Figures  3  and  4  are  represented  the  waiting  time  and  completion 
time  distributions  for  batch  jobs  for  the  same  values  of  N  as  in  Table  I. 

It  appears  clearly  that  increasing  the  number  of  terminals  does  not 
influence  much  the  purely  interactive  components  of  the  system,  but 
disturbs  much  more  the  batch  components,  especially  in  the  critical  region 
N=35  to  47  (note  that  the  system  is  unstable  for  N>48). 

Remark.  To  produce  the  curves  of  Figure  3,  we  have  solved  the  system  (11) 
for  each  value  of  N.  The  system  was  truncated  at  the  percentile  .999  in 
order  to  lose  at  most  a  probability  mass  .001  in  the  tail  of  Wq(x).  We 
have  not  computed  Wq(x)  for  N=47,  because  the  resulting  finite  system  of 
differential  equations  would  have  been  of  order  48*125  which  was  too  large 
for  our  program  to  handle. 

In  Table  II  and  Figures  5  and  6,  we  represent  essentially  the  same 

information  as  in  Table  I  and  Figures  3  and  4  for  different  values  of  w-* 

The  last  row  of  Table  II  indicates  the  maximum  number  N,  of  terminals  that 

max 

may  be  active  without  causing  the  system  to  become  unstable. 

It  appears  that  improving  the  performances  of  the  interactive  components 
of  the  system  has  a  profound  effect  on  the  whole  system,  especially  on  the 
batch  activity.  This  is  reflected  on  each  characteristic,  whether  they  be 
moments  or  percentiles  of  queue  length  or  waiting  time  distributions  (see 


Appendix 


Firstly,  we  prove  that  matrix  H*(z,y,s)  is  a  solution  of  the  following 
matrix  equation: 

(A0+A")H*2+(ArsI-A^-ApH*+zA2+yA^  *  0  (A.l) 

where  A£  is  a  square  matrix  of  order  (N+l)  with  the  unique  non-zero  element 
(^2^01  ”  ^*2* 

Proof.  One  easily  proves  the  following  equations: 

d)t  ho(k,f  ,x)  *  ^  ■j  hQ*hg(  k ,  ,  x ) 

+  NX2h1*hg(k,8.,x) 

+  ul6kl6* 

“  [X  -j+NX^+P  1  ]^q(  k  *  x)  9 

~  x i^i*hg(k » i » x) 

+  (N-l^h^k^x) 

+  M26£l6k0 

“  [x^+( N-l ) *^) * 

^  hj(k,t,x)  =  X]hj*h0(k,e,x) 

+  (N-j)A2hj+1(k,t,x) 
u 2h j_ i  ( k , i-l  ,x ) 

-  [x1  +  (N-j)x2+u2]hj(k,J.,x).  for  2<jsM. 

Remark.  We  define  h^+^(k,t,x)  =  0  for  all  (k,l,x). 

Taking  the  generating  functions  of  the  Laplace-Stieltjes  transforms, 
we  get  (writing  ht  for  ht(z,y,s)): 

J  J 
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x1h*2+NX2h*h*-[\1+NX2+y1]h*-shg+u12  =  0, 

Xih*h*-[X1+(N-l)x2+n2]h^+(H-l)X2h*-sh*+li2y  =  0, 

X 1 h j h0" Cx  t + ( N- j ) x  2+u  2 1 h J+ ( N~ J ) x  2h J+ 1 +u  2h J  _ i - s  h J  =  0. 

This  last  system  becomes  in  matrix  notation: 

A0H*2+A2H*2+AlH*'A2H*'A2H*"sH*+zA2+yA2  =  °- 

One  now  easily  obtained  Formulas  (15)  to  (17)  by  deriving  formula 
(A.l)  respectively  with  respect  to  s,  z,  or  y,  by  right  multiplying  by  e, 
and  by  taking  the  result  for  z=l,  y=l ,  and  s=0. 

The  non-singularity  of  the  matrix  M  follows  from  the  fact  that 
det  M=-v.| (-v2)^(1-p)  1  0  for  p<1. 

In  order  to  obtain  Formulas  (18)  to  (20),  it  is  not  necessary  to 
inverse  the  whole  matrix  M,  it  is  sufficient  to  determine  the  first  two 
rows  and  r^  of  M”1.  They  can  be  obtained  by  solving  two  linear  systems 
very  similar  to  the  system  providing  the  vector  Now  Formulas  (18)  to 
(20)  immeidately  result  from 


m0  = 

mBo =  "-^r 


mBl  =  '-1^-1 » 


nI0  =  ‘^2’ 


nIl  =  ’—1^2* 
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Titles  of  Tables 


I  Characteristics  of  the  system  behaviour  for  different  values 
of  N. 

(X^.030,  X2=.056,  ^  =  .5,  u2=2.4) 

II  Characteristics  of  the  system  behaviour  for  different  values 
of  u2- 

(X^.030,  X2=.056,  y -j  = . 5 »  N=43) 

III  Percentage  in  reduction  of  the  mean  interactive  processing  time, 
the  mean  interactive  sojourn  time,  and  the  mean  batch  sojourn 
time  for  increasing  w2. 
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Figure  1. 
Figure  2. 
Figure  3. 

Figure  4. 

Figure  5. 

Figure  6. 


Titles  of  Figures 


The  model . 

The  traffic  coefficient  p(N)> 

The  waiting  time  distribution  for  batch  jobs.  Different 
values  of  N, 

The  completion  time  distribution  for  batch  jobs.  Different 
values  of  N. 

The  waiting  time  distribution  for  batch  jobs.  Different 
values  of  Ug* 

The  completion  time  distribution  for  batch  jobs.  Different 
values  of 


Table  II 


Characteristics  of  the  system  behaviour 
for  different  values  of 
(A^.030,  A2=.056,  u^.5,  N=43) 


w2 

2.55 

2.50 

2.45 

2.40 

2.35 

2.30 

2.25 

2.20 

2.15 

P 

.401 

.435 

.475 

.522 

.579 

.649 

.734 

.840 

.974 

E[j] 

3.96 

4.21 

4.47 

4.76 

5.08 

5.43 

5.81 

6.23 

6.68 

.5 

3 

3 

4 

4 

4 

5 

5 

6 

6 

.9 

9 

9 

10 

10 

11 

11 

12 

12 

13 

.99 

14 

15 

15 

16 

16 

17 

17 

18 

18 

.999 

18 

19 

19 

20 

20 

21 

21 

22 

22 

E[i  ] 

1.05 

1.23 

1.47 

1.82 

2.33 

3.18 

4.84 

9.38 

69.03 

.5 

1 

1 

1 

1 

1 

2 

3 

6 

48 

.9 

3 

3 

4 

5 

6 

8 

12 

22 

160 

.99 

7 

7 

9 

10 

13 

17 

24 

45 

320 

.999 

10 

11 

13 

16 

19 

25 

37 

68 

480 

E[Wj] 

.5 

1.55 

1.68 

1.83 

1.98 

2.16 

2.36 

2.58 

2.83 

3.11 

1.14 

1.26 

1.40 

1.55 

1.73 

1.93 

2.16 

2.41 

2.70 

.9 

3.74 

3.99 

4.26 

4.55 

4.88 

5.23 

5.61 

6.02 

6.47 

.99 

6.56 

6.91 

7.28 

7.67 

8.10 

8.55 

9.04 

9.56 

10.12 

.999 

8.92 

9.33 

9.77 

10.24 

10.73 

11.26 

11.82 

12.42 

13.06 

E[w] 

21.63 

26.50 

33.26 

43.10 

58.38 

84.48 

136.88 

284.58 

2268.57 

E[wc] 

13.37 

14.50 

15.83 

| 

17.41 

19.31 

21.62 

24.46 

28.00 

32.48 

mlc 

29.00 

31.25 

33.88 

36.98 

40.68 

45.12 

50.54 

57.21 

65.53 

mo 

22.34 

25.66 

30.14 

36.44 

45.89 

61.52 

91.89 

175.10 

1269.67 

ml 

3.98 

4.63 

5.51 

6.75 

8.61 

11.68 

17.66 

34.03 

249.38 

Nmax 

50 

49 

48 

47 

46 

45 

44 

44 

43 

Table  III 


Percentage  in  reduction  of  the  mean  interactive  processing  time 
the  mean  interactive  sojourn  time,  and  the  mean  batch  sojourn 
time,  for  increasing 


mean  int.  proc.  time 
mean  int.  soj.  time 
mean  batch  soj.  time 


-42% 


-95% 


-98.5% 


Figure  1 


N  terminals 


The  model 


The  completion  time  distribution  for  batch  jobs. 
Different  values  of  N. 
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The  completion  time  distribution  for  batch  jobs. 
Different  values  of  u-* 


